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Abstract 

The idea of the universe with a mirror sector having all particles and forces 
identical to those in the familiar sector has been proposed in the context of neutrino 
physics as well as superstring theories. Assuming that all the quark and charged 
lepton masses in the mirror universe are scaled by a common factor, (, as is required 
in one interpretation of the neutrino data, we investigate domains of the parameter 
£ where physical conditions are favorable for cooling in the age of the universe that 
can lead to the formation of compact structures given the initial condition Qg ~ fig 
(B denoting the mirror baryon). In particular we ask whether there is a region in 
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C-space for which primordial Jeans mass mirror clouds cannot cool in the present age 
of the universe. We find that, for most of the area of interest in the parameter space, 
atomic hyperfine structure cooling is effective in a time period short compared to the 
age of the universe but long compared to the free fall time for globular cluster-sized 
objects expected on the basis of simple Jeans length analysis. 
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I. Introduction: 



The idea that our present universe may have a mirror partner evolving with 
identical matter and force content has been invoked for various reasons from time to 
time during the past thirty years. The continuing interest in the phenomenology of 
such models was revived in the late eighties by the observation that the superstring 
theories (Green et al. 1986) lead naturally to such a picture where the known parti- 
cles are accompanied by a duplicate set, with identical properties, but in which the 
two sets have little or no interaction except for that of gravitation (see for exam- 
ple Khlopov et al 1991; Kolb, Turner and Seckel 1985; Hodges 1993 and references 
therein). Most recently this idea has emerged from attempts to understand the 
experimental observations relating to neutrino oscillations (Berezhiani and Moha- 
patra 1995; Foot and Volkas 1995; Berezhiani et al. 1995) such as the solar and 
atmospheric neutrino deficits (Winter 1995) as well as indications from the LSND 
(Athanassopoulos et al. 1995; Hill 1995) data. The idea is particularly attractive in 
any model that requires an ultralight sterile neutrino to fit observations. Our focus 
in this paper will be on the class of models advocated in Berezhiani and Mohapatra 
(1995) and Berezhiani, Mohapatra and Dolgov (1995) where solar neutrino data re- 
quire that the weak scale of the mirror universe be higher than that of the standard 
model. This model (the weak scale asymmetric mirror model, WSAM) has features 
of particular astrophysical interest: no mirror element, except hydrogen, is stable; 
the mirror proton mass is similar to that of the normal proton; and the mirror elec- 
tron masss is perhaps ten to 100 times heavier than that of the normal electron, 
while the value of mirror electric charge is the same as for normal matter. 

In Berezhiani et al. (1995), the possibility was raised that mirror matter might 
be detected by lensing or other gravitational phenomena, based on the mirror struc- 
tures that might form, such as supermassive mirror black holes, mirror machos, etc. 
The purpose of this paper is to study specifically the formation of mirror structures 
in this model so as to permit more detailed investigation of their implications for 
the model. 

The question of initial condensation and star formation of the primordial hy- 
drogen cloud of the familiar universe has been a subject of discussion among cosmol- 
ogists for the obvious reason of understanding the origin of first stars and galaxies 
(see for example, Peebles and Dicke 1968; Silk 1977; Palla et al. 1984; Tegmark et al. 
1996). Its study combines information from atomic physics and statistical physics as 
well as cosmology and chemistry. In our discussion, we take over the relevant ideas 
with appropriately scaled atomic parameters and combine them with the cosmology 
of the hidden mirror world to get a qualitative understanding of structure formation 
in the (exclusively hydrogen) mirror sector. A significant difference is that, whereas 
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dust particles play a major role in understanding most of the detailed structure in 
the familiar sector, the mirror sector for the range of interest of the ( = me/m e 
parameter is devoid of any nuclei heavier than the hydrogen and is thus free of dust 
particles. This to an extent simplifies our discussion and enables us to draw the 
conclusions that we do. Of particular interest to us is the nature and evolution of 
the low mass Population III stars for the visible universe. Our calculations draw 
from the literature on Population III stars, particularly the classic 1968 paper of 
Peebles and Dicke. 

The plan of the paper is as follows. In Section II, we review the necessary ele- 
ments of the weak scale asymmetric mirror model(WSAM). In Section III, we address 
cosmological structure issues, including mirror recombination, primordial hydrogen 
molecule formation and the mirror Jeans mass. It is not surprising that, because of 
the similarity in proton mass, the mirror matter Jeans mass at matter domination 
is about that of a globular cluster; an important difference is that because of the 
greater electron mass, mirror matter recombination precedes both that of normal 
matter and the point of matter domination over radiation. In Section IV we address 
the fate of Jeans mass "mirror globs." We find that for most regions of parameter 
space, the energy loss rate should be sufficient to support considerable condensation. 
We show that, as ( increases, molecular cooling becomes ineffective but cooling by 
the hyperfine structure transition seems to take over sufficiently rapidly as to leave 
in all probability no interesting ( region without a cooling mechanism. In particular, 
we show that the rate of hyperfine cooling is small compared to the free fall rate. In 
Section V we discuss the fate of the globs, including possible scenarios within the 
context of the results of Sections III and IV and possible observational constraints. 

II. Review of the Weak Scale Asymmetric Mirror 
Model 

There exist several experimental indications of non-zero neutrino masses and 
mixings: these include (i) the solar neutrino deficits observed now in four different 
experiments, (ii) deficit of atmospheric neutrinos observed in 1MB, Kamiokande 
and Soudan II experiments and (iii) a preliminary indication of direct accelerator 
neutrino oscillation from the LSND experiment (for a recent review, see Winter 
1995). If all these data are confirmed by the planned experiments such as SNO, 
Borexino and Super Kamiokande as well as the new runs at LSND etc, they will 
have profound implications for the neutrino mass matrices as well as for physics 
beyond the standard model. One simple way to understand all these observations 
is to assume that there exists a fourth ultralight neutrino species beyond the three 
already known (Caldwell and Mohapatra 1993; Peltoniemi and Valle 1993) i.e. z/ e , 
and v T \ in this picture, the solar neutrino puzzle is understood via the oscillations 
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between the v e and the extra neutrino (denoted by u' e ) both of these are assumed 
to have mass in the milli-eV range. The atmospheric neutrino puzzle is solved via 
v ir v t oscillation. The present indications from LSND experiment then imply that 
m v ~ m Vr ~ few eV. Indeed there exist analysis (Primack et al 1994) of the existing 
data on structure in the universe which seem to require a hot dark matter component 
at the level of 20% of the critical density which is fulfilled if the more massive of the 
neutrinos above have a mass of 2.4 eV each. 

The problem with such a picture is that since LEP data implies that only three 
neutrino species can couple to the Z-boson, the fourth neutrino must be singlet (or 
neutral) with respect to the standard model gauge group. Since a standard model 
singlet fermion can in principle acquire an arbitrary mass in the absence of any 
symmetry protecting its lightness, the puzzle arises as to why such a neutrino (i.e. 
u' e ) is so light. It was suggested (Berezhiani and Mohapatra 1995; Foot and Volkas 
1995) that if there is a mirror sector of the universe with identical gauge structure 
and particle content to the standard model but decoupled from the standard model 
particles except thru gravity, then the same symmetry (i.e. B — L) that keeps 
the ordinary neutrinos massless (or light) will have an analog in the mirror sector 
(to be called B' — L') which will keep the v' e also massless or ultralight. It was 
suggested (Berezhiani and Mohapatra 1995) that the lightest neutrino masses arise 
from Planck scale effects, in which case the masses of v e and v' e scale like the square 
of the weak scales in their respective universes (i.e. v^ k and v'^ k ) whereas their 
mixing goes like v' wk /v w k = C From the MSW solution to the solar neutrino puzzle, 
one can then conclude that ( ~ 30. In this picture, the mirror analogs of and v T 
will have masses in the range of few KeV's so that they can constitute warm dark 
matter for the universe. 

It furthermore turns out that since the value of Aq C d in the mirror sector is 
near that in the normal sector whereas the mirror down quark, d! is ( times heavier 
than the ci-quark, we expect m n > — m n ^binding energy of nuclei. As a result, the 
mirror neutron whether bound or free is always unstable (Berezhiani, Dolgov and 
Mohapatra 1995) resulting in the startling conclusion that in the mirror sector the 
only stable atom is the mirror hydrogen atom (H r ). The weak scale asymmetry 
model (WSAM) has also the major difference that the ionization energy of H' is 
perhaps on the order of 400 eV (for ( = 30). 

The basic WSAM picture of the evolution of the mirror sector after the big bang 
has been discussed in Berezhiani, Dolgov and Mohapatra (1995) It was pointed out 
there that consistency with standard big bang nucleosynthesis requires that there be 
asymmetric reheating after inflation in the two universes with T ~ T'/2. In fact, in 
Berezhiani, Dolgov and Mohapatra (1995) we constructed realistic models where the 
asymmetric inflationary reheating is intimately connected with the asymmetry in 
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the electroweak scales. The rest of the history of the mirror sector can be worked out 
until recombination in the mirror sector which takes place when the temperature of 
the universe is around few tens of eV due to higher ionization energy of the mirror 
hydrogen. The nature of structure in the mirror sector depends on the detailed 
dynamics of the mirror sector and is the problem we attempt to tackle in this paper. 

III. The Early Years 

In attempting to study possible structure formation in the mirror universe in 
WSAM, we follow a scenario similar to the familiar sector (Kolb and Turner 1990; 
Weinberg 1972; Peebles 1993 and references therein). For ( mildly greater than 
one (( > 3), mirror recombination takes place before matter dominance. After 
the mirror recombination (at which point the mirror baryon component dominates 
the expression for Jeans length), the mirror universe contains a neutral hydrogen 
cloud with a small fraction of ionized hydrogen and electrons. When the universe 
later becomes matter dominated, the density fluctuations grow with time until the 
expansion of the high density regions stop and violent relaxation takes place leading 
to isolated "globs". The future evolution of the globs is dictated by the rates for 
physical and chemical reactions that determine the extent of energy loss by the 
globs. If the energy loss rate is significant, then condensation to jupiter type objects 
of low mass or supermassive black holes can proceed; otherwise the gas cloud remains 
diffuse making detection difficult. The primary difference between the mirror and 
the familiar sector will be the scale factor ( which will make the reaction rates in 
the mirror universe different. 

We compute: (1) amount of ionization at mirror matter recombination; (2) 
density of primordial mirror hydrogen molecules; (3) Jeans mass and length of a 
mirror matter "glob" at matter domination; and (4) glob parameters after violent 
relaxation. We consider, for definiteness in most places, that the universe has the 
critical density, with normal and mirror baryons each providing 0.05 and mirror 
neutrinos, the model's warm dark matter, the remainder. We take the Hubble 
constant to be 50 km/s-Mpc. We take the ratio of the mirror proton mass to that 
of the normal proton to be a = 1.5 and the ratio of the mirror temperature to the 
normal temperature to be (3 = 1/2. We sometimes keep the ratio of the mirror 
electron mass to that of the normal electron a free parameter (; where a definite 
value is useful, we choose ( = 30. 

In computing recombination in the mirror model, we follow the treatment of 
Kolb and Turner (1990) for normal matter, and reproduce it in order to compare 
the mirror matter calculation with it. We consider the reaction: 

e" H + 7 (1) 
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We define the variable x = B/(2T), where B is 13.6 eV. The reaction rate for (1) is 
then given by: 

M = A x (x^/ml) (2) 



with 



A\ = 47 ^ 2 = 1.46x10 14 cm 3 s 1 MeV' 2 . We have for the equilibrium ionization 



fraction, 

X eq = 0.51^ (^) ^ e~* (3) 

= A 2 x 3 ^e~ x /(n)i 2 h) (4) 

with r] = (tt B h 2 )2.Q8 x 10~ 8 and A 2 = 1.41 x 10 7 . We have for the density of 
electrons, 

n e = ??n 7 X e = ^"^mV^/i 2 ) 1 / 2 (5) 
with A 3 = 2.8 x 10 16 cmT 3 MeV~ 3 . Finally, the time of recombination is 

t R = lH \T /Tf 2 = A 4 x 3 / 2 /(hm 3 ^) (6) 

with A A = 1.55 x 10 10 MeV 3 ' 2 s. 

Setting (crv)n e tR = 1 to find the time at which freezeout occurs gives the 
condition, 

A.A^m'^x-^e-^ 2 = 1 (7) 

1 /2 

Using this one finds, x~ l l A e~ x = 1.56 x 1CT 13 ($f^) • This gives for the fraction f e 
of free electrons 

/e = Xeq = [A^ZaJibH] = 6 X 10 " 5 £^ (2^5) (8) 

This is in rough agreement with Kolb and Turner (1990). 

In the WSAM model, we must modify reaction rates appropriately. Denoting 
all parameters of the mirror sector by a tilde and defining a = mp/m p and (3 = T/T 
and C = — as before, we have for the reaction rate for the analog of reaction (1) 



(av) = AyxW/ml 
7 



(9) 



where we have assumed that the cross section for (1) scales with the square of the 
Bohr radius and x — B/2t. The analogues to Equations (3-6) are 

X eq = A^x^e-t/infh) (10) 



n e = A 3 x- g / 4 m 3 e-*(riy 2 h)/a 1 / 2 (11) 



where A 5 = 9.6 x 10 J MeV 2 s. Setting 1 = (av)n e t^ gives 



™ 1/2 

mat '~ 



X 1/4 e~* = I 1/2 ^ I ( 13 ) 



so that we have 



/ 1/2-1/2 \ 

X^/X R = {A*/ A*) ( (C«/A) ~ 2.5C (14) 

Numerically, Eq(14) corresponds, with our standard parameters (£ = 30, etc), 
to about 10% unrecombined mirror electrons. 

We now use this result to compute primordial formation of mirror matter 
molecules (Hydrogen). Primordial molecule formation has been addressed in detail 
by Lepp and Shull (1984) and in summary by Peebles (1993) H 2 is formed in the 
early universe, before galaxy formation, principally in two catalytic processes 

e~+H^H'+-i (15) 

followed by 

H + H~^H 2 + e- (16) 

and 
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(17) 



followed by 



H + H+ -> H 2 +p (18) 

In each case, the first reaction is slow and the second fast. In each case, the amount 
of H 2 formed is determined by the temperature at which the CBR can no longer 
dissociate the ion (H~ or H 2 )• Lepp and Shull give the two temperatures as 
2.75Kz e ff with z e ff=6A and 190, for H~ and H 2 , respectively; the corresponding 
binding energies are 0.75 eV and 2.65 eV. The amount of H 2 produced, in either 
case, is then 



h = f ( z eff)-^ = ^ v )( z eff)n H (z eff )f e t(z eff ) (19) 

where t is the age of the universe at z e ff. f 2 in Eq (19) rises with z e ff since the z 3 
behavior of nu outweighs the z~ z l 2 behavior of t(z e ff). For familiar H 2 , Lepp and 
Shull show that, although z e ff is larger for Eq (19), the contribution of Eq (15) is 
about three times as big owing (mostly) to the smaller velocity of the proton in Eq 
(17). For mirror matter, however, the velocity of the electron falls like C 2 ? so 
we expect Eq (19) to be a good approximation and dominant for ( > 9. (av) for 
familiar matter, for Eq (19) is a constant (1.4 x 10 _18 cm 3 s _1 is used by Lepp and 
Shull (1984)) at low energy, so a conservative estimate for mirror matter is simply 
to scale it as (~ 2 . z e ff scales as (/{3 since T = f3T, so we take z e ff = 380( . Thus, 
for mirror matter, Eq (17) becomes 

/2 = fe(vv)nH(z eff )t(z eff ) (20) 

with f e = 4 x 10~ 3 C and t(z) = 4 x 10 17 /z 3 / 2 , giving 

f 2 = 1.2 x 10~ 6 C 1/2 (21) 

This is about 60 times the amount for familiar matter, with the same fl, for this pro- 
cess. It would be of interest to compare this rate with results for muonic molecules; 
however the reactions of interest for muon catalysis appear to be only those involv- 
ing atoms, rather than ions, since densities are high (Hughes and Wu 1977). On 
the other hand, the precise numerical result in Eq (19) will not be needed below. In 
Section IV, we will note that the most important factor in molecular cooling, as a 
function of zeta, is that, for given zeta, the temperature be high enough to permit 
exciting molecules out of the rotational ground state. 

Third, we turn to the Jeans mass. We consider, as noted, a three component 
system - massive mirror neutrinos dominating the mass density, together with both 
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normal and mirror baryons. The "baseline value" for matter domination is z=5800; 
and the eigenvalue Jeans equation is (Kolb and Turner, 1990) 

5iv1k 2 - AnGN^jPjdj = (22) 

where Si are the density fluctuations in the various components, (i = v,B,B) pi 
are the mean densities, Vi are the sound velocities in the different media and k is 
the wave number for the density fluctuation. The various velocities are given by 
v p = y/3f/m 9 , = v^T^m^f;^) 1 / 2 , v* B = l/y/Z 
The solution to Eq (22) is 

S t = A/v 2 ; k 2 = 4nGZ jPj /v 2 = (23) 

We note that the other eigenvalues of Eq. (22) are all zero, for any number of 
components, and that the amplitudes in the solution of Eq. (23) are driven by the 
inverse of the sound speeds; the densities do not enter. The evolution of galaxies 
presumably proceeds at leisure as in warm mixed dark matter models; see, for 
example, Colombi, Dodelson and Widrow (1996) and references therein for recent 
detailed discussion. It is clear that, as the coolest component by far, the mirror 
baryons dominate the solution. We find a Jeans length and mass of 

Aj ~ 7 x 10 18 {z 6 /zmd) 1/2 &-j 1/2 cm (24) 

Mj = 2.9 x 10 38 nj 1/2 (z MD /z 6 ) 3 / 2 (25) 

where Zq = 5.6 x 10 3 and Qj = £7^/0.05 The mass in Eq. (25) is solely that of the 
mirror baryons within the Jeans length. The basic result is, not surprisingly, on the 
order of a globular cluster. In Table 1, we collect some of the features of the early 
thermal history. 

Finally, we need to make assumptions about the rate of growth of the inhomo- 
geneity. We adopt a simple parametrization. We assume that the mirror hydrogen 
globular cluster size glob continues to expand with the rest of the universe from the 
time at z=5800 of initiation of matter domination until 



Zstop = 58002m 
10 



(26) 



At that point "violent relaxation" begins. It leads to glob of half the original (z stop ) 
size with temperature determined by the virial theorem. The Mass, radius, Temper- 
ature, and mirror particle density in the glob, and the age of the universe are given, 
in terms of the parameter zm, in Table 2. There will, of course, be a distribution 
in mass of the structures; we restrict this first look at the structure question to just 
those of approximately the Jeans mass. 

The prediction of the mirror matter model is then formation of structures 
roughly on the order of those shown in Table 2. There is a strong dependence on 
the choice of the parameter zm in the table. We, of course, do not know the rate at 
which inhomogeneities in mirror matter will grow, but we do know that they should 
become nonlinear earlier than those of normal matter, since mirror recombination is 
earlier (and indeed before the time of matter domination) and since 5i in Eq (23) is 
significantly larger because of the smaller sound velocity. In the estimates below, we 
shall consider 0.1 and 0.01 values that can't be ruled out for zm, but not necessarily 
maximum and minimum values. 

IV. Energy Loss 

In this section we estimate rates at which the mirror matter globs can lose 
energy by different processes so as to determine how compact present day objects 
might be in order to estimate their susceptibility to gravitaional detection. We 
consider three processes: (1) free- free radiation - as modified by continuing recom- 
bination; (2) the atomic mirror hydrogen hyperfine structure line (or 7 mm line); 
and (3) the ortho-para mirror molecular rotational transition line. The basic tension 
in cooling is that the temperature scale is set by the time of matter domination and 
so is somewhat similar to that for normal matter; but the processes are all strongly 
dependent on the mirror electron mass which is greater by 10 to 100 than for nor- 
mal matter. There is also strong dependence on the parameter zm which, as seen in 
table 2 enters in most quantities of interest. We make estimates mostly on the basis 
of constant (average) densities, but consider where necessary modifications for an 
isothermal sphere. Note that, since in the mirror sector there is no stellar burning, 
there is in all likelihood no reionization and we do not consider such a possibility. 

We begin with the question of glob recombination which will determine the 
availability of mirror electrons for free- free radiation. Using Eq (9), with T as in 
Table 2, we can write for the time rate of change of the mirror electron fraction in 
the glob 



^ = -n H f e (ov) = -f e /t q (27) 
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where t a ~ 10 10 z M 1 sec. The solution is 



fe = fefl/[l + tfe,o/t q ] 



(28) 



We thus see that free electrons will be reduced by at least a factor of ten in time to > 
the current age of the universe, so long as zm is greater than 0.01, and considerably 
more for higher values. We note that the gammas from glob recombination do 
not constitute a thermal energy mechanism for the glob since the Thomson cross 
section is too small for them to heat the remaining electrons, so that they either 
remove binding energy from the glob in the case of recombination to excited states 
or reionize mirror atoms in the case of gammas from recombination to the IS state 
too far from the glob surface. 

The rate at which the electron gas can absorb energy T from an atom is given 



where we have estimated the cross section as being just the mirror bohr area, Eq. 
(27) assumes that there is no problem with radiating the energy; it just computes 
the time needed to transfer the energy from an atom to the free electrons (in fact 
the bremsstrahlung rate will fall as ( increases). The result is that, although fewer 
scatterings are needed as zeta increases, the cross section falls so that the time 
required is nearly independent of zeta. Note that the factor of (m e -/m^) 2 represents 
the random walk in energy space needed. The actual rate of energy loss will be 
still less because of opacity, but Eq (29) shows, taking into account the decrease in 
fe of Eq (28) as z M increases from 0.01, that significant condensation via free-free 
radiation must take more than the age of the universe. Note that, since we have a ~ 
C~ 2 > fe ~ C> v ~ C" 1 ^ 2 i n Eq.(29), this conclusion is independent of (. This argument 
also indicates that Compton cooling (scattering of electrons off CBR photons) (see 
Silk, 1976) cannot cool either in the parameter range under consideration. 

Next, we turn to a particularly interesting possibility, open to mirror hydrogen 
because of the increased energy per transition and the increased transition rate that 
come with the greater electron mass- the atomic hyperfine transition. The time 
needed to radiate an atom's total kinetic energy in units of the hyperfine transition 
energy is given by 



by 




(29) 




(30) 
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where the factor in the first set of parentheses is the number of transitions each 
atom needs to make, while e^ F and t^ F are the energy and inverse radiation rate 
for the case of normal matter. Here we have used the fact that the hyperfine splitting 
goes like m\/m^. We can derive a second characteristic time in addition to that of 
Eq.(30) which is the total time required for the atoms to make the transitions that 
radiate away their kinetic energy. The second characteristic time is that required to 
excite the atoms to the higher energy state so that they can do the radiating. The 
reason for treating the two times separately is that each must be less than the age 
of the universe in order for the glob to cool by this process. The total excitation 
time is given by 



In Eq. (31), we have taken the cross-section to be six Bohr areas (Allison and Del- 
garno 1969). We see from Eq (30) that energy loss from the glob, in less than the age 
of the universe, tends to be possible for small values of z M and large values of ( that 
is, when the kinetic energy to be dissipated is small and the energy dissipated per 
transition is large. We note, however, that our calculation is a simple-minded one; 
we have not considered such possible additional effects as collisional de-excitation. 
We note that for ( > 1000, the energy difference between the two levels would in- 
crease to the eV range thereby making it difficult to excite the atom to the higher 
level (i.e. making t2,E larger), thus disfavoring this as a cooling mechanism. 

Finally, we turn to cooling by molecular radiation. This is the chief source of 
cooling in the Population III globular cluster work, for example of Peebles and Dicke 
(1968) and Tegmark et al. (1996). In the present case its effectiveness is enhanced 
(for ( > 1) by the increased energy carried away in each transition, but is decreased 
by the difficulty of exciting rotational states, since even the lowest-energy excited 
state soon becomes more than the temperature that corresponds to z M = 0.1. We 
can write for the time needed to remove energy 3T/2 from every atom 



where we use again the bohr area for cr, use fi as computed in Section II, and 
have 0.015eVX 2 for AE using the experimental value for the ortho-para transition 
for ( = 1. The lowest excitation energy is three times the size of the "average" 
rotational energy ((m e /mH)E e i); characteristic energy levels are given in Table 3. 




(31) 




(32) 
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The factor A 2 takes account of the fraction of collisions that are sufficiently energetic 
to excite the J=0 to J=l, para to ortho, hydrogen molecule transition, and is given 
by 



/d 3 t;id 3 t;2e-( 2m "i +m ^)/ 2r g(m/3(i;i - v 2 f - AE) 
2 ~ / d 3 v 1 d 3 v 2 e-( 2mv i +mv %V 2T 

Eq. (33) is written for collisions of H atoms with H 2 molecules (hence the factor 
of 2 in the exponent); however A 2 is independent of the masses of the particles and 
depends only on AE/T. Evaluating Eq.(33) numerically, we find, approximately, 
that the time in Eq.(32) is greater than the age of the universe for the curve ts = to 
in Fig.l. That is, the glob can lose energy by molecular radiation for larger values 
of zm and smaller values of ( , since in this case it is the ability to excite that is 
important. The range of ( for which molecular radiation can be effective will be 
further limited by the decreasing (with () cross-section for collision induced radiation 
but it will not be necessary to compute this additional effect nor do we compute the 
extra time needed to excite to J = 2 state so that radiative quadrupole radiation 
can occur without collision. 

We now consider the question of opacity. For the case of the rotational line, 
following Lang (1980) and Rohlfs (1986), we have for the opacity. 



c 2 



K,, — 



8tTZ/ 2 

where A and are given by 

32tt 4 



f 2 n H A[l - e' vlT \4>{v) ~ 3 x KT 18 rW4f (34) 



A = ^_au 3 a%/C 2 (35) 



and 

c / M \ 1 ' 2 



*w = i;(2m) ^ 104/2l/v (36) 



From Eq (34) one can check that 90% of the mass is within an optical depth of 300 
from the surface so that total radiation time does not become a limiting factor in 
energy loss. One can also see from Eq (35) that, for the case of the magnetic dipole 
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hyperfine structure line the opacity will be reduced by at least a factor of 10 5 so 
that the optical depth will be less than one. 

We have found that: (i) free- free radiation cannot give enough energy loss 
to permit collapse for any (, for reasonable zm ', (ii) hyperfine structure radiation 
cannot give such energy loss if ( is too low, for reasonable zm because the rate at 
which energy is radiated is too small and it cannot give such energy loss if zeta is too 
high (over 1000), for reasonable z M) because the rate at which the upper hyperfine 
state is excited becomes too small; and, (iii) finally, molecular radiation cannot give 
such energy loss if zeta is too high (over 5), because the ortho to para molecular 
transition cannot be excited with the available thermal kinetic energy. These results 
are summarized in Figure 1 where we show the three curves 

t 2 ,R = l0 21 z M /( 3 = t =Ax 10 17 sec (37) 



h,E = lO^VVAi = t (38) 

where Ai is the analogue of Eq.(33) for HH scattering. Finally, the condition t 3 = to 
leads to 

z% 2 h. 2 = 1.2 x 10~ 8 (39) 

In Fig 1, the curve left-most at the bottom (t 3 = t ) is that above which the 
glob can cool by means of the ortho-para Hi transition; below it, the temperature of 
the glob after violent relaxation is too low to excite the transition. The second curve 
{t2,R — to) is that below which the glob can radiate away its energy by means of 
the hyperfine transition; above it, too many transitions are needed to radiate away 
the thermal energy in time to- The rightmost curve is that above which the glob 
can excite atoms to the upper hyperfine state a sufficient number of times to radiate 
away the thermal energy; below it the density and excitation rate are too low. In the 
shaded area between the first two curves and under the third curve cooling cannot 
take place. Our result is that cooling is likely possible for most of the region ( < 400 
and zm > 0.01 with the lower bound on logzM rising approximately like 1.6log(. 

V. Where Are They Now? 

Given the above results, we can now ask the likely fates of mirror globs and the 
possibilities for detecting signs of their existence. There are three broad possibilities: 
(i) Puffy Globs; this would be the case if we are in a region of the zm — C plane in 
which the globs cannot lose energy within t . (ii) Cluster Globs; this would be the 
case for a region of the Zm — C plane in which the globs can lose energy sufficiently 
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rapidly, and there is fragmentation to low mass bodies, which would essentially be 
white dwarfs or " small" black holes, since no nuclear burning is possible, (iii) Black 
Globs; this would be the case for a region of the zm — C plane in which there is 
energy loss and no fragmentation, and hence collapse to a supermassive black hole. 

Case (i) will obtain for the shaded regions of Figure 1, providing there are 
not important energy loss mechanisms we have neglected. Case (ii), in principle, 
could obtain for portions of the Zm — C plane in which the energy loss is through 
the hyperfine line. This is because the glob will be optically thin to the hyperfine 
line and hence, if the energy loss were rapid enough, particle velocities will decrease 
and regions of smaller Jeans length will appear. See, for example, the discussion of 
Bohm-Vitense (1989), or the original paper of Hayashi (1966). However, the rate of 
energy loss is never rapid enough from the hyperfine line. That is, we can take the 
total time for energy loss by adding Eqs. (37) and (38), divide by the free fall time, 
(Gp) -1 / 2 , and minimize the ratio with respect to (. We find 

(t 2R +t 2E )/t ff > 10 (40) 

Were the left hand side of Eq. 40 much less than one, we could expect frag- 
mentation through decreasing Jeans mass leading to subregion collapse. Case (ii) 
will also obtain if the mechanism of Peebles and Dicke (1968) for sequential central 
star formation in globular clusters (followed by glob reheating and radiation from 
central objects formed) is applicable, or if other energy loss mechanisms turn out 
to be important. The mechanism of Peebles and Dicke might maintain case (i) for 
the current age of the universe. Knowledge of star formation, for ordinary matter, 
especially formation of Population III stars, is not sufficiently advanced, as far as we 
can determine, to draw any confident conclusions. Case (iii) will, of course, obtain 
where there is energy loss without fragmentation to case (ii). 

Given these three cases, we can look a little further into the fate of the globs 
by asking about glob-glob scattering under the assumption that the globs are bound 
to galaxies. We take a simple model of our galaxy as a sphere of radius 70 kpc, with 
a virial velocity on the order of v g = 10~ 3 c The number of collisions a glob is likely 
to make is then given by 

N c = (N g /V g )a g v g to « 10- 3 ^ 2 (41) 

where we have assumed N g = 10 6 globs in the galaxy, and taken as the cross section 
the area corresponding to the glob radius in Table 2. We have neglected gravitational 
focusing. We can, in addition, ask about the angle of scattering. We have (see, for 
example, Binney and Tremaine 1989) 

tanO « Av/v = 2GM/(bv 2 ) « 0.1z M (42) 
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The results of Eqs. (39) and (40) are that we are near a borderline; we might 
have a lot of scattering, including some large angle scattering that could liberate 
the glob from the galaxy. It would appear unlikely that we would have less than 
on the order of a percent of the globs liberated. Thus, one might consider searches 
for shadow matter globs both inside the galaxy and outside. Searches inside benefit 
from the possible existence of a large number of globs close at hand. Searches outside 
benefit from the large ratio of globs to galaxies, even with only 1% liberation. 

We note that we would expect, in case (ii) above that some percentage of the 
cluster elements would be liberated in glob-glob scattering, but we would consider 
it unlikely that a large fraction of mirror matter goes into making MACHOs. We 
should also note that cases (i) - (iii) above are not necessarily mutually exclusive: 
while there will only be one value of zeta, we expect some range of values of % 
for globs. If near a border, these could include more than one case in the above 
taxonomy. Finally glob-glob mergers could play an important role in determining 
some mirror sector structures and should be considered in more detail. 

Finally, we turn to the question of potential observational constraints. The one 
that comes most immediately to mind is that of lensing, galactic and extragalactic; 
detailed analysis is needed and possible based on the above results, but is beyond 
the scope of the present paper. There are, in addition, other possibilities. One 
interesting way to limit the number of globs that could be resident in the galaxy is 
through the study of wide angle binary star systems. These would be disrupted by 
close enough passage of massive mirror globs in the same way they are disrupted 
by passage of giant molecular clouds as discussed by Bahcall, Hut and Tremaine 
(1984) There is similar application by Hut and Tremaine (1985) to the question 
of loss of the sun's Oort cloud, the source of (long-period) comets. Applying the 
results of Bahcall, Hut and Tremaine, we find a time for (tidal) disruption of wide 
angle binary systems by mirror globs somewhat less than the age of the galaxy. 
However, as emphasized by Weinberg (1985) in addition to the uncertainties in our 
model, there is the problem of not knowing the initial numbers of wide angle binaries 
with enough confidence to say how many have been disrupted. We therefore do not 
consider this a reliable test of the model. We note in passing that there is a bound 
(Anderson et al. 1995) on how much nonluminous matter, including mirror matter 
is allowed to exist around the sun. 

Another fruitful area for detecting mirror globs would be work on limiting 
the number of halo black holes the galaxy can have, including observational work 
such as searching for X-rays from infalling matter. With 10 6 globs in the galaxy, 
there should be considerable contribution to the X-ray background. In a separate, 
straightforward direction, the above rough model of the galaxy with a million globs 
in a 70 kpc sphere would have more than one glob within 1000 pc of the sun. One 
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might be able to look at peculiar motions of stars within such a distance to see 
whether they are consistent with a large, but unseen, nearby mass. 

There would be other gravitational effects on galaxies from mirror globs. Rix 
and Lake (1993), consider limits on 1O 6 M objects. They show that such objects 
could contribute to dark matter for the Milky Way but that 1O 3 M is the maximum 
allowed by relaxational disk heating in two spiral dwarf galaxies. They point out 
that this result would not hold for objects that could be dissipated in the higher rate 
of collision for the dwarf case. Their result argues against black compact globs but 
not against the other cases. It probably does not rule out case (hi), since glob-glob 
scattering might in the dense environment of the dwarf galaxy, evade the constraint 
by the liberation from the galaxy of globs and by evaporation from the globs of some 
of their mirror hydrogen before collapse. 

VI. Conclusion 

We have analysed the cooling mechanisms for a Jeans mass glob in the mirror 
sector as a function of the two free parameters ( and % which respectively denote 
the mass ratio of electrons in two sectors and the ratio of the redshifts at which the 
mirror globs undergo violent relaxation to that of matter domination. We found 
that for 5 < ( < 400, cooling is likely to occur by means of the atomic hyperfine 
structure line. This overlaps with the range for ( (between 10 to 100) required to 
understand the results of the solar neutrino observation. While we show that the 
dominant cooling mechanism in the ( region of interest is too slow to give rise to 
glob fragmentation by itself, we cannot rule out with confidence any of the three 
end-point cases (at the current age of the universe) of: (i) puffy globs not greatly 
condensed, (ii) cluster glob with fragmentation to low mass bodies or (iii) black glob 
with collapse to a single black hole. We noted several possible ways to test and/or 
constrain the mirror idea and the three endpoints observationally including lensing, 
wide angle binaries, effect on relatively local stellar phase space distribution and 
black glob effects such as x-ray production. 
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Table 1 



T 


z 


tu 


n 7 


n B cm 3 


T = 2.7K 





t = 4 x 10 17 sec. 


422cm" 3 


9.4 x 10~ 8 


T R = 0.26eV 


1100 


1.1 x 10 13 sec. 


5.6 x 10 11 


1.3 x 10 2 


T EQ = 1.4 eV 


5.8 x 10 3 


7.0 x 10 11 


8.2 x 10 13 


1.8 x 10 4 


T^=14(C/30) eV 


5.8 x 10 4 (C/30) 


2.9 x 10 10 (30/C) 2 


8.2 x 10 16 (C/30) 3 


1.8 x 10 7 (C/30) 3 



Caption: Early history of mirror matter; here t v stands for the age of the universe. 



Table 2 



M g 


2.9 x 10 38 


Rg 


3.5 x 10 18 /^ M 


T 


1.8z M eV 


Pg 


4 x lO" 19 ^ 


n g 


1.6 x 10 5 z 3 M 


tg 


6 x !0 12 z M 3/2 



Caption: Glob parameters after violent relaxation. 

Table 3 



Energy levels 


Electronic 


vibrational 


Rotational 


Mirror matter energies 


10C eV 


.2C 3/2 eV 


5 x 10~ 3 C 2 eV 



Caption: Mirror matter energy levels 

Figure Caption: The lines correspond to the three cooling conditions t 2t n = to, 
ts = to and ti,E = to, where to is the age of the universe. The shaded areas denote 
the range of the parameters zm and ( where cooling cannot take place before the 
present age of the universe t . 
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